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We show that voids adjacent to radiation damage cascades can be moved in their entirety by 
several lattice spacings. This is done using molecular dynamics cascade simulations in iron at 
energies of 1-5 keV. The effect of this mechanism is studied further using an OKMC code and 
shows enhancement of void diffusion by 2 orders of magnitude from 1x10" cm /s to 3x10"^" cmVs. 
Repeating the study on He bubbles shows that the movement is damped by the presence of helium 
in the void. 

PACS numbers: 61.82.Bg, 61.72.Qq, 61.80.Az 



I. INTRODUCTION 

Steels under high energy neutron bombardment such 
as experienced in fission and fusion reactors are known 
to develop voids and bubbles over time. These have se- 
riously deleterious effects, and a clear understanding of 
how these defects affect the properties of a material is 
crucial to designing better steels. Molecular dynamics 
can be used to gain insight into the atomistic level inter- 
actions between radiation cascades and defects such as 
voids and bubbles. The long term effects of the atomic 
mechanisms can be observed using object kinetic Monte 
Carlo (OKMC) and such a multi-scale approach is be- 
coming increasingly popular. 

A study of very low energy cascades of just a few eV 
impacting on voids has been performed by DubinkJ^. It 
showed that very low doses can lead to the dissolution of 
voids by focusons striking them. This is consistent with 
experimental observation^ that the voids can be reduced 
in size by irradiation under favourable conditions. 

At higher energy PrP et al. studied stability of helium 
bubbles struck by radiation cascades. They find that 
the stability of the bubble depends on its He/ vac ratio. 
Clusters with significantly more He than vacancies tend 
to absorb vacancies from the cascade. Conversely clusters 
with a low density of He tend to lose vacancies during the 
cascade, approaching the stable ratio of near 1:1. 

Further properties of He in iron during cascades (but 
not in the presejjce of voids or bubbles) have been ob- 
tained by Lucas^ et al. They find that the presence of 
substitutional He tends to reduce the number of Frenkel 
pairs produced. Interstitial He increases lattice stress and 
favours the formation of self interstitial atoms (SIA) as 
well as stabilising SIA clusters. They observe increased 
diffusion of He at temperature, which explains how bub- 
bles can form rapidly. 

An interesting mechanism for interstitial cluster for- 
mation during cascades in pure a-iron has been observed 
by CaldePet al. They find that some cascades emit par- 
ticles moving faster then the cascade wave front, which 
cause secondary cascades ahead of the main one. Colli- 
sion of the high density primary wave with the core of 



the secondary cascades can lead to nucleation sites for 
interstitial clusters. This is due to atoms being forcefully 
pushed into the secondary cascade's low density core re- 
gion. 

Under certain irradiation conditions voids can form 
latticed. This effect has been reproduced by 
Heinisch^ et al. in KMC simulations. Their study uses 
clusters of crowdions formed by irradiation and moving 
mostly one dimensionally along close-packed planes to ex- 
plain the alignment. Voids off the lattice are exposed to a 
greater flux of crowdions and so tend to shrink, whereas 
voids on-lattice are shielding each other along the close 
packed planes. The simulations do produce a void lattice 
successfully although the lattice does not show the same 
refinement as experimentally observed. 

This current study looks at the interaction of voids and 
bubbles with cascades and how this affects their centre of 
mass. A mechanism is observed of atomic injection into 
the void by the density wave-front of the cascade. This 
leads to a new void being formed at the cascade core, 
because the atoms absorbed by the original void are not 
available for recombination. The mechanism is explored 
atomistically in molecular dynamics and its consequence 
for diffusion rates is explored via using OKMC. 

II. METHOD 

The molecular dynamics study used a version of 
MOLD^ modified for variable time step cascade simula- 
tions as detailed previouslj^. The Ackland-Mendele\ff^ 
embedded atom potential was used for Fe-Fe interactions. 
The Juslin pair potentiaP^' was used for Fe-He and the 
Becli!^ pair potential for He-He. The Beck He-He po- 
tential was splined to the ZBlP universal potential at 
short range. Helium is a closed-shell atom, and we do 
not expect any chemical bonding, so the pair potential 
formalism is sensible. 

The study consists of triggering radiation cascades in 
a-iron near to a void or bubble. To create stable voids 
or bubbles in a BCC iron lattice, a single atom was re- 
moved near the centre of the system. The system was 
then relaxed and the atom with the highest potential 
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energy removed. Iterating this process produces a stable 
void with faceted sides ahgned to stable lattice directions. 
Stability was tested by running the sample at 900 K for 
20 ps and observing no changes. 

Bubbles were produced by placing helium atoms onto 
the vacant sites in the void. The system was then 
quenched using low temperature MD with a high thermo- 
stat value to limit the velocity of the helium atoms. Once 
the helium has reached a stable configuration the system 
may be thermally equilibrated with standard thermostat- 
ing. 

Cascades were initiated by giving a substantial kinetic 
energy to one iron atom near to the defect cluster, this 
atom is called the primary knock-on atom (pka). Cas- 
cade energies of 1, 2.5 and 5 keV were considered. For 1 
and 2.5 keV cascades a box dimension of 40 (128000 
atoms) was used, with 65 Qq (549250 atoms) for 5 keV. 
We use the NVT canonical ensemble (constant number, 
volume, cell shape, temperature and zero total momen- 
tum) with periodic boundaries. 

The final atomic configurations were analysed by com- 
parison to a perfect lattice aligned to the defective 
lattice'^. The Wigner Seitz cell around each perfect lattice 
point is checked for atoms, allowing easy identification of 
interstitials, vacancies and substitutional defects. De- 
fects were considered to be clustered if they were within 
2"^^ nearest neighbour distance of each other. 

The data was further processed by comparing the near- 
est neighbours of every atom before and after the colli- 
sion, marking atoms with at least 2 new neighbours. This 
identifies the area that has been affected by the cascade, 
as well as being sensitive to correlated motion of atoms. 
Additionally it provides an effective method for visualis- 
ing the cascade by joining a line from the initial to final 
location of each atom (c.f. BacorPet al.). The region cov- 
ered by these lines is the part of the simulation "melted" 
by the cascade and the lines tend to form connected paths 
due to replacement-collision events. By taking the list of 
atoms that have changed neighbour, and filtering it so 
that only members that have moved at least one lattice 
site remain, an estimate of cascade volume is obtained. 
The cascade radius is calculated from this volume. 

The point around which the cascade is centred is iden- 
tified by the centre of mass of those listed atoms with 
at least four neighbours also in the list. This eliminates 
small sub-clusters and replacement-collision chains that 
would otherwise skew the estimate. 



III. RESULTS 

Simulations of a cascade interacting with a void or bub- 
ble in a periodic cubic supercell of a-Fe were performed 
under various conditions. Most pkas were given initial 
velocities close to (111) direction towards the void, with 
the slight deviation to avoid perfect replacement-collision 
chains. This high symmetry direction was chosen in or- 
der to localise the cascad^^^ and give easily-comparable 
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Figure 1: Graph showing movement of void, Ad cluster, (filled 
squares) and bubble (empty triangles) vs. distance to PKA, 
d, for a 30 defect cluster hit at 1 keV. It can be seen that the 
void can move significantly towards the PKA site for certain 
cascade distances. The presence of helium in the void greatly 
damps this efi^ect. 




d (angstroms) 

Figure 2: Graph showing movement of void. Ad cluster, (filled 
squares) and bubble (empty triangles) vs. distance to PKA, d, 
for a 100 defect cluster hit at 5 keV. It can be seen that the 
void can move significantly towards the PKA site for certain 
cascade distances. The presence of helium in the void greatly 
damps this effect. The effect is much more pronounced than 
at 1 keV. 

data across different simulations. 

To explore if the effects observed are dependent on 
the crystal alignment some cascades were also per- 
formed in the (311) direction and in random directions. 
Simulations were performed at 1, 2.5 and 5 keV on 
voids/bubbles of size 30 and 100, with the bubbles at 
ratio 1:1 for vacancy :helium. One set of data is also 
performed at 300 K to explore temperature effects (the 
various simulations along with their damage yields are 
summarised the appendix tables). 

For 300 K simulations little change was seen in the re- 
sults apart from a slight reduction of Frenkel pairs, and so 
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Figure 3: (Color online) The mechanism for cascade-induced 
void movement is shown above for a 5 keV collision with a 
100 vacancy void. The blue (dark grey) dot represents the 
starting positions of the cascades, orange dots (grey) are va- 
cancies, light blue dots (white) are interstitials or filled vacan- 
cies. Small black dots represent the atoms that have moved 
in the cascade. Non-defective atoms are omitted and the full 
box size is not shown. Top left: Void and pka positions and 
PKA direction shown after 2 fs. Top right; Wave of displaced 
atoms starts to form and expand at 76 fs. Bottom left: Full 
extent of wave, encompasses the void, injecting atoms into its 
core at 260 fs. Bottom right: Wave collapses back, but atoms 
are trapped in void, leaving a depletion at cascade core by 
around 1 ps (image taken at 5 ps for clarity). Note if an 
initially-vacant site in the void is filled, it is coloured light 
blue (white). This was done to show the filling of the void. 




Figure 4: (Color online) a) (left): Cascade radius overlaps 
half the void leading to partial filling and splitting of the void. 
Half the void is unaffected, and now forms a smaller void. A 
second void is formed at the depleted cascade core, as several 
atoms from the cascade are now trapped in the near side of 
the old void, preventing their return. The shaded region rep- 
resents atoms that have been displaced by the cascade, and is 
related to the cascade radius. Orange dots (grey) are vacan- 
cies, light blue dots (white) are interstitials or filled vacancies, 
b) (right): Typical final state of a 5 keV cascade impacting 
100 He bubble in a 100 vacancy void. He remains clustered 
in the void. He atoms are white with a dot. 



further temperature studies were not pursued, although 
higher temperatures may have more effect. Results along 
the (311) direction were found to be consistent vi/ith sim- 
ilar studies of the (111) direction. The radii and dam- 
age of the cascades was also consistent across randomly 
launched trajectories. 

Normally, voids move by capture and emission of point 
defect, but interestingly some of our simulations show 
that the void is moved by several angstroms, far more 
than absorption of the defects created by the cascade 
would allow. The process appears to have an upper and 
lower interaction radius, close to the size of the cascade, 
as shown in Fig. [T] & [2j Detailed inspection reveals the 
diffusion process. 

A typical cascade displaces atoms around itself creat- 
ing a density "wave" moving spherically outwards, with a 
depleted region at its core. Normally in a perfect lattice 
this "wave" will collapse back as the energy dissipates, 
causing recombination of vacancies and interstitials at 
the core (Figure [3|. The defect yield (dpa) from a cas- 
cade normally considers only point defects surviving this 
recombination. 

However, when a pre-existing void is adjacent to the 
cascade a more complex mechanism may take place. The 
density "wave" injects atoms into the pre-existing void, 
where they find sites to occupy and so do not return as 
the wave collapses back. This results in the original void 
being filled, and a new void of equal size appearing at the 
cascade core. The outcome of this process is equivalent 
to the diffusion of the entire void by several lattice spac- 
ings. This behaviour was reproducible for voids under 
all conditions considered, providing the cascade overlaps 
the void. 

Figures [T]& [2] show how the distance moved by the cen- 
tre of the void depends on the separation between void 
and cascade. The cascade will typically begin a few lat- 
tice spacings from the starting position of the pka. At 
very short distances the PKA can pass directly through 
the void, which is why there is a minimum interaction 
threshold for PKA distance. Data suggests interaction 
thresholds of 12 A and 25 A with peak interactions 
around 18 A for 1 keV (Fig. [I]). For 2.5 keV the inter- 
action starts at 18 A , lasting till 40 A and has peak at 
31 A. For 5 keV the interaction starts around 12 A ends 
at 38 A and had peak interaction at 31 A (Fig. [2|. 

For IkeV cascades on bubbles of size 30 almost no mo- 
tion of the cluster is observed even for high cascade over- 
lap. The slight movement seen (Fig. [T]) is due to the side 
of the bubble nearest the cascade absorbing a few vacan- 
cies, shifting the centre of the cluster. On average the 
bubble grew by one vacancy and almost all the helium 
remained clustered in it. 

For a void of size 100 impacted at 5 keV, significant 
movement of the void towards the cascade was observed 
(Fig. [2| . Void size tended to decrease slightly on average 
to 84 vacancies. This change is largely due to the void 
being split, as shown in Fig. [4^, which occurs in 3 simu- 
lations. Ignoring the split samples void size tends to be 
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around 88 after cascade interaction. 

Similar results to the 1 keV bubble case were observed 
for a 5 keV cascade impacting a helium bubble of size 
100. Fig. |4]d shows a typical bubble after impact. The 
helium reduces the mobility of the cluster (Fig. [2]) and the 
helium remains clustered in the bubble. Cascade radius is 
unaffected, although the number of Frenkel pairs formed 
seems to increase slightly (table Occasional shrinkage 
by up to 5 vacancies is observed at longer distance due 
to interstitial absorption from the cascade periphery. 



Table I: The Frenkel pair count, Nf, and cascade radius, 
for various simulations. E is the pka energy. 



E (keV) 


rc (A) 


type 


Nf 


1 


7.3 


void 


5.5 


1 


7.2 


bubble 


6.2 


2.5 


10.8 


void 


11.2 


5 


15.1 


void 


16.7 


5 


15.3 


bubble 


19.1 



Table II: Model parameters: Table of migration energy, Em, 
and dissociation energy of point defects from clusters, Ed- 
Formation energies: EJ = 2.07 and Ej = 3.77 eV. Et is the 
binding energy of a cluster with two defects, E^ = 0.30 and 

,,2/3_(„_i)2/3 



El = 0.80 eV. All diffusion is 3D. 



22/3-1 



Type 


E^ (cV) 


Ed (cV) 


V 


0.67 






0.62 


0.97 


Vs 


0.35 


1.04 


Vi 


0.48 


1.29 


V„,n>4: 


immobile 


E^:^' + Ej + « - Ej)g{n) 


I 


0.34 


h 


0.42 


1.14 


h 


0.43 


1.26 


h 


0.43 


1.26 


I„,n> 4 


immobile 


E:^^ + Ej + {Et - E'f)g{n) 



IV. OKMC RATES 

To explore the significance of this mechanism in terms 
of enhancing void diffusion, an object kinetic Monte 
Carlo (OKMC) model was written. It implements the 
standard resident-time algorithm, also called the kinetic 
Monte Carlo or Gillespie algorithrrP^. The established 
mechanism for void diffusion is by point defect absorption 
and emission events (later referred to as standard diffu- 
sion) and so our new mechanism's diffusion is compared 
to this process. The OKMC model includes vacancy clus- 
ters, interstitial clusters, defect creation from cascade de- 
bris and direct cascade-void interactions. It allows diffu- 
sion of small clusters and thermal dissociation of point 



defects from clusters following the energy barriers given 
by the Fu^^ et al. study as used in Caturla's^'' et al. 
model (see table [n]). The capture radii of clusters follows 
the model of Caturla^® et al. 



z 



I/V 



where ?'/„,v„ are the capture radius for interstitial and 
vacancy clusters respectively. Zj = 1.15 is the bias fac- 
tor accounting for the increased strain field observed for 
interstitials. Vl is the atomic volume and rg = 3.3 A. 

Open absorbing boundaries were used for this study, 
simulating a single void at the centre of a single grain. 
Radiation cascades were introduced into the sample with 
the maximum interaction radius taken from the data sets 
in the current study (Figsjl] & [2]). The cascade debris 
is calculated by the Bacorfi^t al. empirical formula for 
iron: Nf = 5.67E'^ ™ , where Nf is the number of Frenkel 
pairs and E is the damage energy in keV. Cascades were 
implemented in a shell like manner, with vacancies ran- 
domly placed in a smaller radius, and surrounded by a 
shell of interstitials. Cluster interactions were carefully 
performed to take into account the correct change in clus- 
ter position for absorption and emission events, as this 
will affect the diffusion. 

If the cascade is sufhciently close to the void, the void 
is shifted to the centre of the cascade region as seen in 
the MD simulations. By checking the interaction volume 
of the cascade for the largest void within distance and 
updating according to: 



if '^cascade ~f '^void '^i 



case— vol 



d then TyQid — "^cascade 



f cascade IS the cascade radius, r^oi^; the void radius, 
is the distance from the cascade centre to the 



^case—vote 



void centre, and r„ 



i/r. 



cascade 



is the position vector for 



the void/cascade. This mechanism was disabled to sim- 
ulate the standard diffusion. 

The trend of both diffusion mechanisms with tempera- 
ture, dose rate, cascade energy and system size was stud- 
ied. The contribution from each process was calculated 
by separating out the jumps due to each mechanism. 
Studies were performed with default values: Tempera- 
ture of 500 K, dpa/s of lxlO-^ cubic cell of 100 A and 
cascade energy of 10 keV on a void of 100 vacancies. 
Only one parameter was varied from default at a time, 
to gain information on system trends without performing 
an exhaustive search of the parameter space. Parameter 
ranges considered: Temperature of 400-600 K, dpa/s of 
10-8-10-2, cubic cell of 50-200 A and cascade energy of 
1-20 keV. Under the default parameters the standard dif- 
fusion was found to be IxlO-^^ cm^/s while the diffusion 
due to the new mechanism is SxlO-^'^cm^/s , two orders 
of magnitude larger. 

For fixed dpa rate, the standard diffusion rate is found 
to not be strongly dependent on cascade energy. The 
mechanism only depends on the number of free point de- 
fects and increasing cascade energy decreases the num- 
ber of cascades to conserve dpa rate, however, the new 
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mechanism is highly dependent on cascade energy, be- 
cause this sets the interaction radius. 

System size does not strongly affect the new mecha- 
nism. It does affect the standard mechanism as the only 
initial sinks in the system are the void and cell walls. 
Increasing cell size gives defects more chance to be ab- 
sorbed by the void, instead of the walls and so effectively 
increases the void capture radius. This leads to a slight 
increase in diffusivity. 

Rates for both mechanisms increase approximately lin- 
early with dose rate so this parameter is not directly crit- 
ical, although it is important in setting the relative time 
scale of the system as compared to dissociation events. 

Neither mechanism is strongly temperature dependent. 
The standard mechanism lacks temperature dependency 
because the diffusion of point defects is much faster than 
the time between cascades. Temperature will only be- 
come critical at lower T when diffusion time of point 
defects is comparable to the time between cascades. 

The new mechanism is found to be dominant over the 
standard mechanism at all parameters considered apart 
from for energies of 2 keV or below, when the interaction 
radius becomes too small to have significant effect. Cas- 
cade energy is the most significant factor, with diffusion 
of the new mechanism rising to 1x10"^® cm^/s at 20 keV. 



V. CONCLUSIONS 

We have shown that the migration rates of voids in ion- 
irradiated environments arc much higher than expected 
from simple point defect migrations. Cascades occurring 
adjacent to voids may cause a jump of the entire void by 
several lattice spacings, if the cascade is within appro- 
priate interaction thresholds. This motion is due to the 
injection of atoms by the cascade into the original void, 
trapping them and leaving a new void at the cascade 
core. The effect becomes more significant as cascade en- 
ergy increases (even at fixed dpa) because it depends on 
cascade radius. This jump mechanism has been shown 
in OKMC simulations to be a more significant mecha- 
nism than diffusion due to absorption and emission of 
point-defects. 

The mechanism is suppressed when helium is present 
in the void. 
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Appendix A: Simulation data 



Table III: Supplementary table of simulations: Voids. Nf is the number of EYenkel pairs formed. Vyoid is the number of 
vacancies in the largest cluster. Tc is the radius of the cascade. 



E (keV) 


setup 


simulations defect-size 


Nf 


^void 


re (A) 


1 


(111) 


17 


30 


5.3 


28.8 


7.1 


1 


random 


48 


30 


5.8 


29.8 


7.2 


1 


(311) 


10 


30 


5.9 


29.6 


7.5 


1 


(111), T=300 K 


14 


30 


4.6 


28.1 


7.6 


2.5 


(111) 


12 


30 


14.2 


30.8 


11.0 


2.5 


random 


16 


30 


9.9 


26.8 


10.6 


2.5 


(111) 


9 


100 


9.4 


89.6 


10.9 


5 


(111) 


12 


30 


18.8 


27.9 


15.2 


5 


(111) 


16 


100 


14.6 


84.3 


15.1 



Table IV: Supplementary table of simulations: Bubbles. Nf is the number of Prenkel pairs formed. Hebub is the helium 
remaining in the bubble after the cascade. Hesub is the helium that has become substitutional in the lattice. Hei„t is the 
helium that has become interstitial in the lattice, rc is the radius of the cascade. Vi,„i, is the number of vacancies in the bubble. 
Hsbub is the number of helium atoms in the bubble. 



E (keV) 


set up 


simulations He:vac bubble 


Nf 




H&bub 




HCint 


rc (A) 


1 


(111) 


14 29:30 


6.5 


31.1 


28.3 


0.6 


0.1 


7.4 


1 


random 


46 29:30 


6.1 


30.2 


28.9 


0.1 





7.1 


5 


(111) 


12 100:100 


19.1 


102.2 


100 








15.3 



r 



